TUHEP-2000-01 
PDK-744 

Search for Nucleon Decay with Final States 

£^Tf\ urf, and Vtt^-''^ Using Soudan 2 



D. Wall^'", W.W.M. Allison^, G.J. Alner^ D.S. Ayres\ W.L. Barrett^ 

C. Bode^, P.M. Border^, J.H. Cobb^, R. Cotton^ H. Courant^, J. Chung^ 

D. M. Demuth^, T.H. Fields\ H.R. Gallagher^, M.C. Goodman^, R. Gran^, 
O : T. Joffe-MinorS T. Kafl^a^ S.M.S. Kasahara^, W. Leeson\ P.J. Litchfield^ 
O : W.A. Mann^ M.L. Marshak^, R.H. Milbum^ W.H. Miller^, L. Mualem^, 
^ A. Napier^ W.P. 01iver^ G.F. Pearce^ E.A. Peterson^, D.A. Petyt^ 

. L.E. Price^, K. Ruddick^, M. Sanchez^, J. Schneps^, A. Stassinakis^, 

^ : J.L. Throng S.P. Wakely^, N. West^ 

in ■ 

^ ! ^Argonne National Laboratory, Argonne, IL 60439 

>■ I "^University of Minnesota, Minneapolis, MN 55455 

I ^Department of Physics, University of Oxford, Oxford 0X1 3RH, UK 
' '^Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 

g '■ ^ Tufts University, Medford, MA 02155 

! ^ Western Washington University, Bellingham, WA 98225 

o ; 

■ "^Now at Sapient Corporation, One Memorial Drive, Cambridge, MA 02142 

1) ; 

(D , 

^ . Abstract 

' We have searched for nucleon decay into five two-body final states using 

a 4.4 kiloton-year fiducial exposure of the Soudan 2 iron tracking calorimeter. 
For proton decay into the fully visible final states fi'^-q^ and we observe 

zero and one event, respectively, that satisfy our search criteria for nucleon 
decay. The lifetime lower limits {t/B) thus implied are 89 x 10'^'^ years and 
81 X lO^'' years at 90% confidence level. For neutron decay into urf , we obtain 
the lifetime lower limit 71 x lO^'^ years. Limits are also reported for neutron 
decay into Ftt^, and for proton decay into Ftt"'". 



PACS numbers: 11.30.Fs, 12.20.Fv, ll.30.Pb, 14.20.Dh 



I. INTRODUCTION 



A. SUSY predictions for nucleon decay 



The decay of the nucleon is a possible experimental window into fundamental processes 
at high mass scales. Calculations carried out in supersymmetric (SUSY) grand unified 
theories (GUTs), such as SUSY SU(5) and SUSY SO(IO), have become increasingly more 
detailed, and predicted lifetimes are within a range which may be accessible to experiment. 
In contrast to their non-SUSY precursors, SUSY GUT models generally favor decay modes 
with final state K° or K"^ mesons. Some of these models also predict that other decay modes 
may exhibit a significant branching fraction. The latter modes include the lepton plus eta 
meson and the uir modes. In SUSY SU(5) models, lifetimes predicted for i'^t] and Vt] are 
longer than those for z/K or Vir modes by factors which vary from 2-3 to one or two orders 
of magnitude . Eta mode lifetimes which are hundreds of times longer than z/K or i'^K 
lifetimes have also been predicted in SUSY SO(IO) [^0. However, it has been proposed 
that in SUSY GUTs such as SO(IO), there exists a new set of color triplets and thereby 
a new source of d=5 operators which may allow p — > i~^7] to become prominent (along with 
VK~^ and Vn^). In SUSY GUTs the Vtt modes are usually predicted to have decay lifetimes 
which are longer than VK modes. In specific models, however, it is possible for interference 
between third and second generation amphtudes to alter this situation so that Vti modes 
become dominant 

We have previously searched for nucleon decay final states with strangeness using the 
Soudan 2 iron tracking calorimeter and reported the results as lifetime lower limits . In 
this work, we extend our investigations to include the t^rf, vrf ^ Vtt^, and Vtt~^ final states. 



B. Previous searches for £^r], vrp , and vtt modes 



Two-body i^'^r] and I/vr decays of nucleons were regarded as interesting in the context of 
non-SUSY GUT models, and experimental searches were reported at various times during 
the past two decades. Multiple ring topologies were investigated for l^rj, Vi], and z77r° final 
states by the 1MB and Kamiokande water Cherenkov experiments p!0|-pr^. In the HPW 
water Cherenkov experiment, events exhibiting a two muon decay signature were examined 
for compatibility with p —>■ ^^r], with t] decaying to tt+tt^tt^ or 7r^7r^7 ||TB|]. The t^rj, Vrj, 
and Vtt modes were also investigated using the planar iron tracking calorimeter of Frejus 
I^,[TB|. More recently, substantial improvement of lifetime limits for the two-body eta modes 



has been reported, based upon the 7.6 kton-year exposure of the IMB-3 water Cherenkov 
detector |16|. Generally, these experiments observed zero or modest numbers of candidate 
events. Understandably, the more loosely constrained n Vrj and n V-k^ modes were 
found to exhibit more background and consequently to yield less stringent lifetime lower 
limits. 



II. DETECTOR AND EVENT SAMPLES 
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A. Central Detector and Active Shield 



Soudan 2 is a massive 963 (809 fiducial) metric-ton iron tracking calorimeter. It operates 
as a slow-drift, fine-grained, time projection chamber of honeycomb-lattice geometry and 
with dE/dx imaging of tracks and showers. As such, it differs considerably from water 
Cherenkov detectors and from planar iron tracking calorimeters. 

The construction and performance of the Soudan calorimeter are described in previous 
publications [|17|. In brief, charged particles are imaged via one meter long, slightly con- 
ductive, plastic drift tubes. Electrons liberated by throughgoing charged particles drift to 
the tube ends under the action of a voltage gradient applied along the tubes. The tubes 
are sandwiched between mylar sheets so as to comprise a "bandolier" assembly. Corrugated 
steel sheets with interleaved bandolier are then stacked to form a massive lattice. A stack is 
packaged with wireplanes and cathode pickup strips at the drift tube ends and surrounded 
by thin steel skins which provide the gas enclosure. The resulting assemblies are 1 x 1 x 2.5 m, 
4.3 ton calorimeter modules. The tracking calorimeter is constructed building-block fashion, 
with contiguous "walls" which are two modules high and eight modules wide. The calorime- 
ter is surrounded on all sides by a double-layer, cavern-liner proportional tube active shield 
of 1700 m^ total tracking area 0]. Shield augmentations in the form of additional single 
layers and a double-layer top cover of proportional tubes have been operational since 1995. 

For the task of searching for nucleon decay in a variety of final states, the Soudan 2 
calorimeter provides imaging of non-relativistic as well as relativistic tracks, and its sys- 
tematics are different from those of other experiments. The detector is located at a depth 
of 713 meters (2070 mwe) in the Tower-Soudan Underground Mine State Park in northern 
Minnesota. Data-taking is still underway, having commenced in April 1989 when the total 
mass was 275 tons, and continued as more modules were installed. Operation with the 
calorimeter at full mass was first achieved in November 1993. The analysis reported in this 
work is based on a total (fiducial) exposure of 5.52 (4.41) kiloton-years, obtained with data 
taken through October 1998. 



B. Data and Monte Carlo event samples 

Four distinct event samples have been assembled for this analysis; two of the samples 
are generated as Monte Carlo (MC) simulations with full detector response, and two are 
comprised of data events recorded in the experiment. These samples, details of which have 
been presented elsewhere are: 

i) Nucleon decay MO. For each of the nucleon decay processes investigated, we 
use MC simulations which track all final state particles through the detector 
geometry. Electronic hits are generated, and detector noise is included, in the 
same format as with real data. 

a) Atmospheric neutrino MC: Events are generated representing charged current 
and neutral current reactions which are initiated by the flux of atmospheric 
neutrinos. The neutrino MC program is based upon the flux calculation of the 
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Bartol group for the Soudan site |T^. The neutrino MC sample analyzed here 
corresponds to an exposure of 24.0 fiducial kton-years. 



Hi) Rock data: We analyze events for which the veto shield recorded coincident, 
double-layer hits. Such events originate with inelastic cosmic-ray muon interac- 
tions in the cavern rock surrounding the detector. These shield-tagged "rock" 
events provide a reference sample by which to gauge cosmic-ray induced back- 
ground events which are included in the "gold data". The latter events arise 
either from shield inefficiency or from instances where an energetic neutral par- 
ticle emerged from the cavern walls with no accompanying charged particles. 

iv) Gold data: Data events for which the cavern-liner active shield array was quiet 
during the allowed time window comprise our "gold event" sample. These events 
are mostly reactions initiated by atmospheric neutrinos but may contain nucleon 
decays, as well as the muon induced rock events with no shield hits (described 
above). The events of interest for our study are those with a "multiprong" 
topology, having two or more particles (other than recoil nucleons) emerging from 
the primary vertex. They are distinct from the more populous single-track and 
single-shower events which originate predominantly from neutrino quasi-elastic 
reactions. 

Events of all four samples are required to be fully contained in a fiducial volume which 
is everywhere 20 cm or more from the outer surfaces of the calorimeter. In order for an 
event to be included in any of the above four samples it must survive the selections of a 
standard processing chain. At the head of the chain are the requirements imposed by the 
hardware trigger; events satisfying these requirements are subjected to a containment filter 
code. Events which survive are subjected to two successive scanning passes carried out by 
physicists; each pass involves three complete, independent scans. In the first scan pass, 
events with multiprong topologies used in this study are found with an overall efficiency 
of 0.98 

-OM- Descriptions of our hardware and software selections, and of our scanning 
procedures are given elsewhere 0,^,^. Our procedures ensure that MC simulation events 
pass through identical or otherwise equivalent (e.g. the hardware trigger is implemented via 
software for MC events) steps in the chain ||^. 

Events which emerge from the second scan pass with a multiprong topology assignment 
are then reconstructed using an interactive graphics package. The set of reconstructed tracks 
and showers which comprises each event is subsequently entered into an event summary file 
from which kinematic quantities are calculated. 

C. Background from rock events and from neutrino interactions 

Among the fully contained, multiprong events of our gold data sample, a small contribu- 
tion may arise from cosmic ray induced rock events. These are events initiated by neutrons 
emerging from the cavern rock which impinge upon the central detector and which are un- 
accompanied by coincident hits in the surrounding shield array. In contrast to neutrino 
interactions or to nucleon decay, these shield-quiet rock events tend to occur at relatively 
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shallow penetration depths into the calorimeter. Their depth distribution, and their dis- 
tribution in visible energy and in other variables, can be inferred from rock events tagged 
by coincident shield hits. To estimate the amount of zero- shield-hit rock background in 
our gold multiprong sample, event distributions of gold data have been fitted to neutrino 
MC and shield-tagged rock samples using a multivariate discriminant analysis P,|21|. We 



find that, the fraction / of all multiprong rock events which have at least one shield hit is 
/ = 0.94 ±0.04. 

For each individual nucleon decay channel, the same event selections applied to data 
multiprongs are also applied to the shield-tagged rock multiprongs. The zero-shield-hit 
rock background is then estimated as the product (1 — /)// = 0.064 times the number of 
shield-tagged rock events which satisfy the channel selections. 

To calculate rates for background events in our nucleon decay search which arise from 
interactions of atmospheric neutrinos in the detector, we use our realistic neutrino MC 
simulation which is based upon atmospheric fluxes with null oscillations. The neutrino 
Monte Carlo program has been described in previous publications 0^ • 



During the past decade, evidence for depletion of the atmospheric muon-neutrino flux as 
described by z/^ —>■ oscillations has become increasingly extensive; especially compelling 
are observations by Super-Kamiokande of zenith angle distortions in fluxes of both sub- 



GeV and multi-GeV muon neutrinos The disappearance of t'^ flavor neutrinos by 

oscillations effectively reduces background in our search arising from {v^ + charged- 
current reactions, and so a correction for this effect to our null oscillation estimates is 
warranted. To implement a correction we assume, as indicated by recent data, that is an 
active neutrino which is not z/g (i.e. Vx = Vt) @- Then, atmospheric neutrino oscillations 
do not affect background arising from {u^ + Ve) charged-current reactions, nor do they 
affect background from neutral current reactions (initiated by any flavor). Consequently, to 
correct for z/^ flavor disappearance, we simply scale the number of {v^ + V^) charged-current 
background events estimated from the null oscillation MC for each nucleon decay channel 
by the v^/ve flavor ratio R measured in the Soudan 2 experiment: R = 0.64 ± 0.13 |2^. As 
noted below in Sections IV and V, and in Table II, this correction yields small reductions in 
null oscillation background rates for p fi~^ri and p Vn^ final states. For p — > e'^r] and 
n — s> Vn^ channels however, the neutrino oscillation correction has negligible effect. 



III. NUCLEON DECAY SIMULATION AND SEARCH METHOD 

A. Event generation 

For each nucleon decay final state, a Monte Carlo sample is created and processed as 



described in Section |II B| . This sample is then used to determine the topological and kine- 
matic properties that differentiate it from the atmospheric neutrino and the rock event 
backgrounds. 

For each final state, about 500 events are generated and embedded into pulser trigger 
events from the detector which are taken at regular intervals throughout the exposure. In 
this way, the detector's evolving size and the background from natural radioactivity and 
cosmic rays are accurately incorporated into the simulation. 
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For two-body decay of a nucleoli at rest, the final state momenta are uniquely determined. 
However, Fermi motion within parent iron and other nuclei of the calorimeter medium smears 
these momenta and thereby complicates final state identification. In our simulations, Fermi 
motion effects are modelled according to the parameterization of Ref. p5| . 

A Monte Carlo event for p fir], t] —>■ 77, which illustrates the search topology for this 
mode, is shown in Fig. |l|. Here, the /i"^ appears as a single non-scattering track. The two 
gammas from rj decay give rise to two showers which are spatially well separated and which 
point to the event vertex which is also the origin of the muon. The fi~^ endpoint decay - 
discernible in approximately 60% of events - appears in Fig. |I] as extra ionization "hits" 
(tube crossings) in the vicinity of the track's range-out point. 

In each nucleon decay mode, the meson can undergo intranuclear rescattering within 
the parent nucleus. For nuclei which have interior as well as surface nucleons {A > 12), 
there is significant probability for event final states to be altered. For final state pions (in 
neutrino MC interactions as well as in nucleon decay), intranuclear rescattering is treated 



using a phenomenological cascade model ||26|. Parameters of the model were set by requiring 
that the threshold pion production observed in z/^-deuteron {A = 2) and z/^-neon {A = 20) 
reactions be reproduced |^. To account for the intranuclear rescattering of rj mesons 
within iron nuclei we have adopted the survival fraction of 0.57 estimated by the Frejus 
collaboration |n| using a detailed balance calculation. This value for the 77° survival is 
similar to the survival fraction of 0.52 which we calculate for the vr" of n ^ z77r° in iron 
nuclei. We note that a recent analysis of 77 meson photoproduction in nuclei finds the 
in-medium r^N cross-section to be nearly independent of 77 momentum between and 500 
MeV/c 0]. 



B. Search contour in the Minv versus |jOriet| plane 

Two quantities that are useful for selection of nucleon decay candidates and the rejection 
of backgrounds are the invariant mass, M^v, and the magnitude of the net three-momentum, 
\pnet\, of the reconstructed final state. 

We create a scatter plot of invariant mass versus net event momentum for the recon- 
structed final states for each simulation. Then we choose a region in this plane whose 
boundary defines a kinematical selection which can be applied to the data and to the back- 
ground samples. We observe that, for most nucleon decay final states, event distributions 
in each of these variables are approximately Gaussian. Consequently, the density distri- 
bution of points on the invariant mass versus momentum plane can be well represented 
by a bi-variate Gaussian probability distribution function. A detailed description of this 
construction is given elsewhere P,p9|. 

Projections of bi-variate Gaussian surfaces onto the Minv versus Ipnetl plane for the eta 
modes of this search are shown in Figs. 2-7. In these figures, the distribution of a nucleon 
decay event population is depicted using five nested, elliptical boundaries. Proceeding out- 
ward from the innermost contour, the bounded regions contain respectively 10%, 30%, 50%, 
70%, and 90% of the MC nucleon decay sample. From the five regions delineated we choose 
the 90%-of-sample contour - the outermost, solid-curve ellipse in Figs. 2-7 - to define our 
"primary" kinematic selection. That is, a candidate nucleon decay event has reconstructed 
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{Minv, \Pnet\) values which lie within the outer contour. The interior contours (dashed el- 
hpses in Figs. 2-7) are helpful to gauge whether an event sample as a whole exhibits the 
kinematics of nucleon decay, and for this reason we display them. For our search involving 
the ?7 — >■ 77 decay, containment within the 90% contour is our sole kinematic constraint; no 
subsidiary condition, such as a cut on the 77 invariant mass, has been used. 

All nucleon decay simulation events are subjected to the triggering requirements, the 
detector containment requirements, and to the scanning and topology cuts. Subsequently, 
kinematic cuts and additional topology cuts, designed to reduce the background on a mode- 
by-mode basis, are applied. The cumulative effect of these selections is to reduce the overall 
detection efficiency significantly. Typically, the product of the triggering, containment, and 
scanning selections reduces the survival fraction to below 30% for any mode. Table | shows 
the successive survival fractions (including 77 decay branching ratios) for all modes studied. 
Note that the effects of intranuclear rescattering processes (INR) are included in the survival 
fractions listed in Table |. For comparison, the rightmost column of Table | shows our 
estimates of ex BR in the absence of INR within the iron and lighter nuclei (atomic masses 
12 to 56) of the calorimeter medium. 

Sections |I^ and below describe the analysis for each nucleon decay mode that we 
studied. In each case, the particular characteristics of the decay are reviewed and the 
kinematic cuts designed to eliminate background are presented. The signal and background 
events which pass the cuts are then tallied. Then, using the mode detection efficiency and 
detector total exposure, a lifetime limit r/B at 90% confidence level is calculated. 
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TABLE I. Survival fractions through event selections applied in succession, for the MC generated 
nucleon decay final states of this study. The column labeled BR refers to the 77 or tt decay. The effect of 
intranuclear rescattering (INR) can be seen by comparing the two rightmost columns. 



IV. SEARCH FOR NUCLEON DECAY INTO ir] AND T7ri 
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A. Search for p — > t^rj 

We have searched for proton decay into ix^rj and e^r] and for neutron decay into Vrj. 
The decay sequences involving the two largest branching modes of the r/, namely rj —>■ 'yy 
and 1] 7r°7r'^7r°, have been investigated. The results from both of these r] decay modes are 
included in the calculation of the limits for all three nucleon decay modes. We also explored 
the possibility of inclusion of r] —>■ 7r"'"7r~7r° and rj — > 7r+7r~7 into our search, but we found 
the resulting events to be difficult to identify with topology criteria. Since the potential gain 
is modest, we have not included these processes. 

To calculate partial lifetime lower limits, t/B, we use a formalism common to previous 
analyses by us and by the Frejus tracking calorimeter experiment [Q. The mode p 



— ^ y^^?7, the first to be discussed below and which involves two daughter processes for the rj, 
provides an example whose generalization is straightforward: 

^)p^M+'?o > J^p^ J^f ^ \ ■ (4.1) 

Here N^f^n) = 2.87(3.15) x 10'^^ protons (neutrons) in a kiloton of the Soudan 2 detector, 
Tf = 5.52 kiloton years is the full detector exposure, and x Bi{ri) are the selection 
efficiencies given in Table |. The Hi are the constrained 90% CL upper limits on the numbers 
of observed events, and are found by solving the equation 

^ . . E:r=o ^n7=o Pjni, h + fii)Pin,, h + /i2) .... 

subject to the constraint 

ei X Bi{r]) ^ 62 X B2{v) ^ ELi Q x Bjjr]) 

In Eq. ( [4.2|) , P{n,fi) is the Poisson function, e~'^fi^/n\, and the bi are the estimated back- 
grounds. 



We have searched for proton decay into fi~^ri, where the t] decays into two photons or 
into three tt*^ mesons. With final states involving either of these t] decays, the /i"*" is easily 
distinguished from the rj decay products in the Soudan 2 tracking calorimeter (see Fig. p. 
For the p — *■ lu.'^r] mode, the /i"*" momentum is 296 MeV/c in the proton rest frame. In the 
laboratory frame, /i"*" momenta are smeared about the two-body nominal value as a result 
of the Fermi motion of bound protons. 

For p —>■ fi~^i], 1] 77, the kinematic region in the invariant mass versus net momentum 
plane which contains 90% of the reconstructed MC proton decay events, is delineated by 
the outer contour displayed in Fig. |^. Since the entire final state is visible, the invariant 
mass will approximate the nucleon mass and the net event momentum will be distributed 
in accordance with the convolution of Fermi motion with detector resolution. Fig. ^ shows 
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MC Soudan2 Data 
Run 33907 Event 1048 
07-Feb-1992 16:16:42.86 



93l 



L - 1. 



117 z (cm) 209 

FIG. 1. Monte Carlo event with full detector response for proton decay p ji^rp^if 77. Here, 
the anode [X) versus drift time {Z) projection has been selected from the three scanning views with an- 
ode-time, cathode-time, and anode-cathode projections. The open square denoted "MARK 1" indicates the 
reconstructed primary vertex. 
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FIG. 2. For p l^'^Tj^ V ~* 77i event samples and the kinematic selection contour in the Minv versus 
\pnet\ plane. The distributions show single track plus two shower events (solid circles) for a) the proton 
decay simulation, b) the atmospheric neutrino Monte Carlo, c) the rock events, and d) the gold data. 
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FIG. 3. For p M^*?; ^ ~^ 7r"7r°7r°, event samples and kinematic selection contour in the Mi„„ versus 
|p„et| plane. The distributions show one track plus 3-6 shower events (solid circles) for a) the proton decay 
simulation, b) the atmospheric neutrino MC, c) the rock events, and d) the gold data events. 
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the distribution of events within this plane, for the proton decay simulation sample, for 
the atmospheric neutrino MC sample, and for the rock and gold data samples. The total 
background which arises from neutrino interactions and from the cosmic ray induced rock 
events is 0.9 + 0.1 = 1.0 events. The neutrino background originates from assorted multiple 
pion production channels which, after correction for z/^ oscillations, is comprised of equal 
portions of z/^ and z/g inelastic charged- current events. 

No events are observed to pass the primary kinematic contour. Due to the small back- 
ground expectation, no additional kinematic constraints (e.g. on the momentum of the rj or 
the fi'^) were applied. With an overall detection efficiency of 6.9% we establish a limit of 
t/B > 48 X 10^° years at 90% CL for this decay sequence. 

In the other decay sequence of p f^~^Vy namely r] vr^vr'^vr^, we search for a decay 
signature of one track and three to six decay showers from the three vr'^s. The vertex is 
defined as the end of the track from which the showers emerge. In this case, the rate of the 
rock background is small and the background expectation of 0.5 events arises almost entirely 
from atmospheric neutrino interactions. Fig. ^ shows the kinematic search regions in the 
Minv versus \pnet\ plaue, together with the event distributions from the four samples. The 
detection efficiency for this decay sequence is 5.9%. No data events are observed to pass the 
kinematic cuts. The limit for this decay sequence is calculated to he t/B > 41 x 10^° years 
at 90% CL. Combining the two daughter processes we obtain an overall limit of t/B > 
89 X 10=^° years at 90% CL. 

p — ^ e+ry 

For proton decay into e'^rj with rj — >■ 77, we searched for an event topology of three 
distinct showers emanating from a common vertex. We found it was not advantageous to 
distinguish the primary positron shower from the decay photons of the r]; rather, an overall 
search in the Minv versus \pnet\ plane was conducted. Fig. H shows all three-shower events 
and their relation to the kinematic cut region for the proton decay MC, atmospheric neutrino 
Monte Carlo, rock, and gold data samples. The total background estimate is 0.9 events, of 
which 0.7 events are calculated to be initiated by atmospheric neutrinos. The background 
events are mostly due to z/g and Vf, inelastic charged current interactions. The number of 
nucleon decay candidates observed is one event. From this we deduce a partial lifetime limit 
of r/fi > 38 X 10^0 years at 90% CL. 

For the r] —>■ n^TT^Tc^ mode, the signal topology is taken to be four to seven showers, 
all emerging from a common vertex. The detection efficiency is calculated to be 6.5%. As 
with 7] 77, no attempt is made to identify which of the several showers in the event was 
from the prompt positron. Distributions for this mode are depicted in Fig. |^. Most of the 
background expectation of 0.6 events can be attributed to Ue charged current multiple 7r° 
production. No data events were observed to pass the kinematic cuts for this mode, and we 
obtain a limit of r/i? > 48 x 10^*^ years at 90% CL. Combining the two r] decay modes gives 
an overall limit for p e^r] of r/B > 81 x 10^° years at 90% CL. 
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n — ^ vq 



Neutron decay n — > z/r] with 77 ^ 77 involves a two-shower final state originating from 
a "sharp" invariant mass of 547 MeV/c^. We define for n — Vrj, rj 77, a kinematically- 
allowed elliptical region in the m^^ versus |pnet| plane, as shown in Fig. |[ The four diplots 
of Fig. ^ show the effect of the kinematic contour selection on the neutron decay simulation 
(Fig. ^) and on the two-shower events from the atmospheric neutrino MC, from the rock 
event sample, and from the gold data. The kinematic contour selection is satisfied by 90% of 
the MC Vt] events, whereas the neutrino and rock backgrounds are almost entirely eliminated. 
With the loss due to intranuclear rescattering, our detection efficiency for n Vr], r/ ^ 77 
becomes 7%. For this mode, there are no candidate events. The background from neutrino 
and from rock events is estimated to be 1.7 events. Nearly half of the neutrino-induced 
background arises from neutral-current production of rj and vr^'s. A significant contribution 
also arises from z/g charged current single 7r° production events in which one of the showers 
is not discernible in scanning. For n — > Vr]^^ we obtain a partial lifetime lower limit of 
t/B > 53 X 10^0 years at 90% CL. 

For the other r] decay mode considered, rj n^TT^iT^, the signal topology is four to six 
showers emerging from a common vertex. The detection efficiency for n — >■ T/rj^o^o^o is 5.4% 
including the intranuclear correction. The distributions in the plane of Minv versus \pnet\ are 
depicted in Fig. 0. From Fig. |^, the total neutrino background is estimated to be 1.5 events, 
of which nearly half is Ue charged current interactions with multiple pion production. The 
remaining background is due to inelastic neutral current interactions. Another 0.6 events 
are expected from rock events, bringing the total background expectation to 2.0 events. We 
observe two data events and this gives a partial lifetime limit of r/i? > 22 x 10^° years at 90% 
CL. Combining the two submodes together yields a limit for n ^ Vi] of t / B > 71 x 10^° years 
at 90% CL. 



V. SEARCH FOR N utt^ AND P Vtt+ 

We have searched for the two-body nucleon decay modes which yield a final-state 
(anti-)neutrino together with a vr-meson. For the neutron decay into z/vr", the observable final 
state consists of two photon showers having a restricted invariant mass, and discrimination 
from background can be carried out similarly to n — z/r], 77 77. 

For proton decay into z/tt"*", the observable final state is simply a single charged track 
which can scatter and/or range to stopping with an endpoint decay {tc^ — >■ /x"*" — *■ e"*"). A 
restricted allowed range of pion momentum is implied by the two-body nature of this proton 
decay. However, the background from quasi-elastic z/^ and neutrinos with unobserved 
proton and neutron recoils is substantial, since muon tracks are indistinguishable from vr"*" 
tracks in the Soudan calorimeter, unless of course a distinct scatter is present. Consequently, 
for the p — > z/TT"*" search reported below, we require that an endpoint decay be present on 
candidate single track events. This requirement discriminates against the background from 
quasi-elastic z/^n — > /i"p reactions (though z7^p /i+n background still remains). 
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n — ^ VTv 



Neutron decay into Vtt'^ yields a two-shower final state. The invariant mass values ob- 
tained by reconstructing event shower pairs from the nucleon decay MC simulation comprise 
the ^-coordinates of the events shown in Fig. |^a. The 77 mass distribution peaks close to the 
7r° mass, however it is rather broad (a mean of 142 MeV/c^ with a sigma of 45 MeV/c^ from 
a Gaussian fit). For z/vr^, a bi-variate Gaussian does not provide a good characterization of 
the event distribution in the Minv versus Ipnetl plane. For the purpose of accommodating 
the data, a two-dimensional Gaussian was multiplied by a sigmoid function (a smoothed 
step function). In fitting the data, parameters governing the orientation of the sigmoid and 
the slope of its step were allowed to vary, in addition to the parameters of the Gaussian. 
The result of this fit is the set of contours shown in Fig. |^. As an additional constraint, we 
require that the two showers have an opening angle less than 90°. Events which are thereby 
eliminated are shown as open circles in Fig. ^ 

In our n Vn^ simulation, 86% of events lie within the outer contour of Fig. j^a and also 
satisfy the 77 opening angle requirement. The overall detection efficiency for this mode is 
11%. For a simulation in which the intranuclear effects are not included we find a detection 
efficiency of 21%. The bulk of the INR losses can be attributed to intranuclear absorption 
and inelastic scattering processes which lower the trigger rate for the simulation from 87% 
to 45%. A similar INR effect is observed for p Vn^ . 

For n Vn^ we calculate that 3.8 background events are expected, with 0.9 events 
from rock processes and 2.9 events from neutrino background. Inspection of the neutrino 
background reveals that 75% of the events are neutral-current inelastic single ir^ production 
events in which the recoil baryon is either a neutron that escapes detection or a proton that 
is produced below threshold. The remaining 25% of the background is due to z/g charged- 
current single charged pion production events in which the pion is misidentified as a shower 
and recoil baryons are undetected. In the gold data we observe 4 candidate events. The 
background-subtracted lifetime lower limit at 90% CL is then t/B > 39 x 10^^ years. 



The two-body decay p Vn^, of an unbound stationary proton would produce a 7r+ 
with a momentum of 459 MeV/c. In the Soudan 2 calorimeter, however, protons are mostly 
to be found within iron nuclei. Pions which are created inside such a nucleus can undergo 
intranuclear rescattering before they emerge; and, even if they emerge with their identity 
intact, the 7r+ mesons traverse a dense medium where they can undergo (further) large 
energy degradation due to scattering processes. The net result is that 7r+ mesons of this 
decay mode exhibit only half of their initial momentum on average, as can be seen from the 
result of our full simulation shown in Fig. For our search we require candidate events 
to have one and only one track (no recoil proton or neutron) with ionization compatible 
with a pion or muon mass assignment. The pion momentum as reconstructed from range is 
required to fall within 140 to 420 MeV/c. Additionally, the track is required to have a visible 
endpoint decay consisting of two or more decay shower hits. These cuts reduce the final-state 
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detection efficiency by roughly 50% but remove enough of the neutrino background rate to 
make a search feasible. The total detection efficiency for this mode is 4.6%. 

To gauge the effects on efficiency of intranuclear rescattering within the parent nuclei 
of the calorimeter medium, a proton decay simulation that does not include intranuclear 
effects has been compared with the full simulation. We observe that the trigger efficiency 
increases from 47% to 72% in the absence of intranuclear effects. Additionally, the average 
momentum of the reconstructed pions of the simulation increases from 284 MeV/c to 356 
MeV/ c. Evidently, part of the discrepancy between the predicted two body decay momentum 
of 459 MeV/c and the average momentum of the reconstructed pion tracks of the simulation 
can be attributed to intranuclear scattering. Hadronic scattering processes in the detector 
medium account for the remaining difference. Finally, in the absence of intranuclear effects, 
the overall efficiency for p — z/Tr"*" would be increased from 4.6% to 8.8%. 

With the cuts optimized as described above, we observe 6 candidate events and we 
estimate the background to be 10.5 events in the absence of oscillations by atmospheric 
neutrinos. However, in the presence of atmospheric z/^ v-r oscillations, our background 
estimate must be scaled by the Soudan-2 flavor ratio (see Section IIC) and is thereby reduced 
to 7.7 events. The lifetime lower limit at 90% CL is then r/i? > 16 x 10^" years. 



VI. UNCERTAINTIES AND LIFETIME LIMITS 

The lifetime lower limits reported here are affected by statistical and systematic uncer- 
tainties which arise with detection of each nucleon decay final state and with background 
estimation. The various error sources, and the corresponding fractional variation Atn/tn- 
thereby introduced, are similar to those detailed for our lepton + K° modes search in Soudan 
2 (see Ref. |p. Sect. V). The one exception lies with treatment of intranuclear rescattering 
losses within parent nuclei for the rj and vr modes studied here. In contrast to produced 
K°s (strangeness = +1), eta mesons and pions may have sizeable rescattering probability 
for which there is also significant uncertainty. Based upon uncertainties which arise in our 
phenomenological cascade model [^], with extrapolation to heavier nuclei of pion produc- 



tion observed in z/^ deuteron {A = 2) and neon (A = 20) reactions [^], we estimate an 
uncertainty of 30% for our rescattering treatment. This error is to be added in quadrature 
to the errors (see below) listed for individual channel detection efficiencies, e x BR, in Table 
I. As can be seen from equation (|4.1|) , the INR uncertainty enters directly into At^v/tat via 
the detection efficiencies e^. 

For each channel there is accumulated error on the survival efficiency through selections 
imposed by triggering, software ffitering, scanning, and kinematic cuts; this can be as large 
as 18%, as indicated by the next-to-rightmost column in Table |. For the utt modes, and also 
for Vf], errors enter the lifetime limit through the estimates of background from atmospheric 
neutrino events and from cosmic ray induced rock events. Propagation of background errors 



through relations (|4.1| ), (|4.2| ), and (|4.3| ) for individual channels gives Atn/tn < 20%. We 
conclude that the uncertainty Atat/t/v on the lifetime lower limits reported in this work 
may be as large as 40%. Of course, comparable uncertainties also apply to other published 
limits on the nucleon lifetime. 
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VII. SUMMARY AND CONCLUSIONS 



A search for five distinct lepton + i]{tt) nucleon decay channels has been carried out 
using a fiducial exposure of 4.4 kiloton-years recorded by the Soudan 2 fine-grained iron 
tracking calorimeter. The modes considered are among those proposed in supersymmetric 
grand unified models. For each mode, cuts have been designed which minimize cosmic-ray 
neutrino-, photon- and hadron-induced background while maintaining sufficient detection 
efficiency to allow a sensitive search. From among all of the lepton + pseudoscalar meson 
modes investigated, zero or small numbers of candidate events are observed; in every case 
the occurrence of candidates is compatible with expectations for background. A summary 
of our partial lifetime lower limits t/B at 90% CL obtained with each channel is given in 
Table 0. 
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TABLE XL Background-subtracted lifetime lower limits at 90% confidence level from Soudan 2. Atmo- 
spheric depletion due to neutrino oscillations has discernible effect for neutrino background to p — > I77r+, 
and also to p ^ /x+?7; background rates for null oscillations are shown in parentheses. 



A comparison of our current Soudan 2 limits with results published by the Frejus, 
Kamiokande, and IMB-3 experiments is presented in Table |T|. For the two-body lepton- 
plus-eta decay modes p Z^^??, P — >■ e^?], and n urj, we observe zero, one, and two 
candidates respectively, with comparable expectations for background processes. Our re- 
sulting lifetime lower limits are compatible with those listed by the Particle Data Group ||3^ 
and are the most stringent limits achieved using iron calorimeters. 
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t/B (10'^" years) 



Decay Mode 



Soudan 2 



Frejus [|l4|,|l5| 



Kamiokande [12| 



IMB-3 in 



p — > e'^rj 

n — > Fry 

n — > Fvr'' 

p — > Vtt~^ 



89 
81 
71 
39 

20 



26 
44 
29 
13 
10 



69 
140 
54 
100 

25 



126 
313 
158 
112 
10 



TABLE III. Comparison of 90% confidence level background-subtracted limits in nucleon decay exper- 
iments. 
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